PHASE DIAGRAM OF THE SPIN-1 ANISOTROPIC HEISENBERG MODEL 

WITH A SINGLE-ION ANISOTROPY 

J. Dely, jan.dely@upjs.sk, J. Strecka, jozef.strecka@upjs.sk and L. Canova, lucia.canova@upjs.sk 
Department of Theoretical Physics and Astrophysics, Faculty of Science, P. J. Safarik University, 
Park Angelinum 9, 040 01 Kosice, Slovak Republic 

ABSTRACT. The spin-1 anisotropic Heisenberg model with a single-ion anisotropy is studied using the 
Oguchi's pair approximation. Although the theory is developed for lattices with general coordination number, 
we treat in detail the three-dimensional lattice with the lowest coordination number, i.e. diamond lattice, where 
the critical and tricritical behavior of the system is analyzed as a function of both the single-ion anisotropy and 
exchange anisotropy. 



INTRODUCTION 

The investigation of critical behavior of quantum 
Heisenberg systems, which are useful for the descrip- 
tion of many magnetic materials, belongs among the 
most plentiful areas of statistical physics. In recent 
years, the special attention has been devoted to the 
low dimensional spin-1 Heisenberg models [1-4]. Pan 
and Wang analyzed the effect of a single-ion anisotropy 
on the ground-state properties of the spin-1 Heisen- 
berg model for a general lattice [5]. However, as far as 
we know, similar treatment of thermodynamic proper- 
ties at non-zero temperatures has not been done up to 
now. Therefore, the purpose of this paper is to study 
the criticality and tricriticality of the spin-1 anisotropic 
Heisenberg model with a single-ion anisotropy for lat- 
tices with general coordination number. 

MODEL AND ITS SOLUTION 

The Hamiltonian of the spin-1 Heisenberg system 
in a presence of the single-ion anisotropy D and the 
external magnetic field h is described as 

H = -Y,j[A{s^s^ + sfsy) + s:s;] 

~DY^{Stf^hY,St, (1) 

i i 

where the first term represents the anisotropic ex- 
change interaction. J is the exchange coupling con- 
stant (J > 0) restricted to the nearest-neighbour 
pairs of spins and A is the exchange anisotropy pa- 
rameter (A=0 and A=l correspond to the Ising and 
isotropic Heisenberg models, respectively). Finally, S] 
(7 — X, y, z) are the components of the spin-1 operator 
at sites i. 

Since the inplementation of single-ion anisotropy D 
in the system under investigation can potentially lead 
to first-order phase transitions, we need to know the 
expression of free energy in order to distinguish the 
stable and unstable magnetic phases. Due to this fact 
we adopt the simple Oguchi's pair approximation (OA) 
[B], which is superior to the mean- field method. Thus, 
the Hamiltonian (1) in the OA for a cluster with two 



spins is given by 

= - J[A(5f 5J + 5f 5p + 

-D\{Stf + {S^f]~K}{St + SI), (2) 

in which h^f is the effective field, acting on both the 
spins of cluster, of the form hf,f — J{q — l)m + h. 
q is the coordination number of the lattice and m is 
the magnetization per site, i.e. m = (^(S'f -I- Sj)). 

The Hamiltonian Hfj^ in the matrix representa- 
tion calculated within a standard basis set of func- 
tions |S'f,S'|) {S^ = ±1,0 and = ±1,0) can be 
written in the form of 9 x 9 matrix. Then it is quite 
straightforward to obtain nine eigenvalues by diagonal- 
izing the matrix of H^^. Thus, the partition function 
Z = Tiijexp{—(3H^^) has the following form 

Z = e^(-^+2^) + 2p/('™)cosh(2/3;ie/) 

^2e/3(^-'^/2)cosh[/3y(iJ~V2)^T2J2A2] 
±4e''^cosh(/3JA)cosh(/3/ie/), (3) 

where (3 = l/ksT is the reciprocal temperature {ks is 
the Boltzmann constant). The knowledge of the par- 
tition function allows us to express the relations for 
all thermodymanic quantities. The free energy of the 
two-spin cluster is then defined as 

/ = -/3-ilnZ + J(g- 1)™^, (4) 

wheares the magnetization m per site is obtained by 
minimizing the free energy (4) and is given by 

m = ^{2e'^(^+2-°^sinh(2/3/ie/) 

+2e'^^cosh(/3JA)sinh(/3/ie/)}. (5) 

Through the use of Eqs. (4) and (5) we can completely 
analyze the phase diagrams of the present model. 

Indeed, close to the second-order phase transition 
from the ordered phase (m ^ 0) to the paramagnetic 
one (m = 0) the magnetization m is very small, for 
h — Q, and one may expand Eq. (5) into the form 

m = am + hrTn? + crrt' + (6) 

where the coefficients a, b, and c are functions of 
r, J, D, A, and q. So, hereby the second-order tran- 
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point is determined by conditions a = 1, b = 0, and 
c < 0. On the other hand, the first-order phase tran- 
sitions between different phases can be obtain from a 
comparison of their free energies. 

RESULTS AND DISCUSSION 

Now, let us study the phase diagram of the system 
and the thermal variations of magnetization for the 
case of diamond lattice (g = 4) in zero external mag- 
netic field. The choice of this lattice is motivated by 
the fact that the diamond lattice is a three-dimensional 
lattice with the lowest coordination number and hence, 
one should expect the most obvious impact of quantum 
fluctuations on the cooperative ordering. 
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Fig. 1. The phase diagram of the spin-1 anisotropic 

Heiscnbcrg model. 

The phase diagram in the reduced units d = D/J 
and t = ksT/J for different values of A is shown 
in Fig. 1. The solid and dashed lines represent the 
second- and first-order phase transitions, respectively. 
The black circles denote the positions of tricritical 
points (TCP). As one can sec from this figure, the 
temperature of second-order phase transitions gradu- 
ally declines with decrease of d, for any value of A, 
until the TCP is reached, where the phase transitions 
change from the second-order to first-order ones. By 
further decreasing of d the correspondent first-order 
phase transition lines fall smoothly to zero. It is note- 
worthy that the temperature of second-order phase 
transitions, in the limit d — > oo, does not depend on the 
anisotropy parameter A and is equal to f* = 4.5392. 
Hereafter, it is evident from Fig. 1 that the present sys- 
tem exhibits the highest critical temperatures in the 
Ising limit (A = 0). On the other hand, the exchange 
anisotropy strengthening gradually decreases the criti- 
cal temperature as a result of raising quantum fluctua- 
tions. Finally, an increase of the anisotropy parameter 
A causes a shift of d-coordinate of tricritical point (dt) 
to the higher values. To illustrate the effect of single- 
ion anisotropy on the phase transitions we show in 



the case of isotropic model (A = 1) for some appropri- 
ate values of d. It is obvious from the figure that for 
d = —1.65 and —1.68 the magnetization falls smoothly 
to zero when temperature approaches the critical tem- 
perature tc- This behavior of magnetization is typical 
for the continuous phase transitions and persists un- 
til d > dt {dt = -1.6801 for A = 1). On the other 
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Fig. 2. The temperature behavior of magnetization of the 
spin-1 isotropic Heisenberg model for several values of d. 

The dashed lines typify the discontinuity of magnetization 
at the temperature of first-order phase transition. 

hand, further reducing of parameter d induces the oc- 
curence of discontinuity in the magnetization at the 
critical temperature tc what is characteristic for the 
first-order phase transitions. It is noteworthy that this 
discontinuity increases when the single-ion anisotropy 
d further retreats from the value dt- 

CONCLUSIONS 

In this paper, the phase diagram of the spin-1 

anisotropic Heisenberg model with the uniaxial single- 
ion anisotropy is examined within the framework of the 
Oguchi's pair approximation. We have demonstrated 
that both the exchange anisotropy A as well as single- 
ion anisotropy d have a significant influence on the 
criticality and tricriticality of the model under investi- 
gation. Finally, it should be mentioned that the more 
comprehensive results of our investigation will be pub- 
lished in the near future. 
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